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I. INTRODUCTION. 



This paper is an additional step in our exploration of the gauge invariant variational 
approach suggested in 0. This approach was developed as an attempt to study analyti- 
cally the nonperturbative infrared dynamics of nonabelian gauge theory. The idea is that 
keeping gauge invariance exact at all stages of the calculation may be crucial for consis- 
tent understanding of the ground state structure of QCD. One therefore tries to account 
for the generation of nonperturbative QCD condensates by introducing them as variational 
parameters in the explicitly gauge invariant (Gaussian like) trial wave functional. 

So far this variational approach has been applied to pure gluodynamics. There it seems to 
capture many of the essential features; mass generation, formation of the gluon condensate, 

and asymptotic freedom, ^-fl- Instanton transitions have been analyzed in || and are 
identified with the saddle points in the integration over the gauge group which projects the 
Gaussian wave functional onto the gauge invariant physical Hilbert space. Interestingly we 
found that in the best variational state instantons of large size are suppressed and the large 
size instanton problem that arises in the standard WKB calculation is completely avoided. 
The dynamically generated mass present in the best variational state stabilizes instantons 
at a size p ~ 1/M. 

Although there remain many open questions regarding the reliability of these calcula- 
tions, it is interesting to see how a simple Ansatz of the same type would do in a theory 
with fermions. In particular one would like to see whether including fermions leads within 
the gauge invariant variational framework to the generation of fermionic condensates and 
chiral symmetry breaking. 

As a first step towards this goal, in this paper we will study the simple toy model, 
1+1 dimensional Abelian theory with massless fermions, the so-called Schwinger model ||. 
Our aim is not to further the understanding of the model itself (which has been solved 
exactly many times by different methods and is well understood) but to develop a workable 
generalization of our method to fermionic theory and test it in the simplest possible but 
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nontrivial setting. 

One of the most striking similarities between the Schwinger model and QCD is that in 
neither theory are the asymptotic states represented by the fields of the Lagrangian. The 
Schwinger model, with massless fermions, exhibits complete screening of charges and has 
one neutral asymptotic state; a 'meson' of mass e/y/n. The presence and interconnectedness 
of the axial anomaly, instantons and the massive asymptotic states of the theory has made 
the Schwinger model a rich field of study. 

The outline of this paper is as follows. In Section 2 we generalize the gauge invariant 
Gaussian variational Ansatz to include fermions and perform the variational calculation. We 
find that the best variational state in the fermionic sector does not contain new nontrivial 
structure. In Section 3 we cross check our calculation using the bosonized form of the model. 
We find that our best trial state is in fact the exact ground state of the theory. In Section 
4 we conclude by calculating the fermionic condensate. We recover the correct value of 
the condensate and show thereby that the chiral anomaly is correctly represented in our 
calculation. 

II. VARIATIONAL CALCULATION IN THE SCHWINGER MODEL. 

The Lagrangian density for the Schwinger model reads 



The theory is super-renormalizable and the coupling constant, e, has mass dimension +1. 
The 2x2 Dirac matrices obey the usual anti-commutation relation, {7^,7^} = 2g^ u , where 
spatial indices run from to 1 and g 00 = 1 We write the Schwinger model in the Hamiltonian 
formalism in the temporal gauge, A = 0. This leaves one degree of freedom for the photon 
field and all spatial indices are now suppressed; A±(x) = A(x) and d\ = d. The Hamiltonian 
is, 




(1) 
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with the usual canonical relations, 



[A(x),E(y)]=i6(x-y) 



E(x) = A(x), 



(3) 



[ip{x), ir^(y)] =iS(x- y), ir^x) = iip\x). 

The gauge invariance of a physical state, in this case the vacuum wavefunctional, is ensured 
by requiring that it satisfies Gauss' law, 



As an Ansatz for the vacuum wave functional of the theory we take the product of the 
vacuum wavefunctionals in the gauge and fermion sectors, projected onto the gauge invariant 
subspace. Following |J for the functional in the gauge sector we take a Gaussian with an 
arbitrary width G(p) which is to be treated as a functional variational parameter. 

Next we have to address the Ansatz in the fermion sector. Since we may expect that the 
effects of chiral symmetry breaking show up in the fermion sector as dynamical generation 
of mass, we propose as an Ansatz for the ground state a Dirac vacuum of free massive 
fermions, which we denote by \m). The mass, or rather the mass function m(p), can depend 
on momentum and is to be treated as an additional variational function. In the absence 
of the gauge fields the fermionic operators that annihilate this state are related to the 
original fermionic operators by the well known Bogolyubov transformation. The state itself 
is somewhat similar to the BCS state. The important difference is that the mixing is not 
between particles and antiparticles, but rather between left- and right- moving particles, 
so that the possible condensates in this state preserve the charge symmetry but break the 
chiral symmetry instead. 

The product of the gauge and fermionic vacua is 



The vacuum wavefunctional in the gauge theory must of course be gauge invariant. This is 
achieved by projecting ([5]) onto the gauge invariant Hilbert space. To do this note that (^) 
transforms under the gauge transformation generated by Gauss' law as 




(4) 
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=exp J dx 4>{x)Q{x)^ ^[A,^j] 

= exp^-iy dxdyA' t '(x)G~ 1 (x - y)A 4> {y) S j exp ^-«e ^ dx^ij^fij^ijj |m) . 

where, 

A+(x) = A(x) - d(f)(x) , (6) 

and 4> is the gauge function. The gauge invariant Ansatz for the vacuum wavefunctional is 
therefore written by integrating ^ over all possible gauge transformations, 



V[A,ip] = J D<p exp ^~ J dxdyA 4 '{x)G-\x-y)A 4 '(y) 

X evU[*MH*W*))M. (7) 



In this formalism, one calculates expectation values of local operators with the Ansatz 
for the ground state. Since we are interested in the calculation of physical quantities the 
operators we shall consider are gauge invariant. The expectation of any such operator O 
can be written as 

(0(A,$,^)) = ^ J D(f> DA (m\ exp (-± J dx dy A{x)G~\x - y)A(y^j 
x0(Afi) ex p("^ / dxdyA^{x)G-\x-y)A^y) ) j 
^^ejdxM^H. (8) 
First let us consider the norm of the state, = 1. We first integrate over A, 
Z = Z A j D(f) exp ^~ J dxdyd(j)(x)G- 1 {x-y)d(f)(y)^ 

Z A = Det 1/2 (2G~ 1 ) . (10) 

The integration over the fermions is a little less straightforward. Remembering that \m) is 
the vacuum of a free fermionic Hamiltonian with mass m one can use the standard trick to 



5 



FIG. 1. Shows the quadratic term resulting from the expansion in small e<p. Solid lines rep- 
resent a fermion loop and the dashed lines represent the gauge parameter. Associated with each 
vertex is the factor — icy°5(t). 

express the fermionic matrix element as a path integral. The gauge transformation operator 
in such a representation then has to be understood as an insertion at fixed "time". This 
insertion generates an interaction between the fermions and the gauge parameter, 



/ 



Dip Dip exp 



iSf — ie / dt dx 5{t)p{x)ip{x)'~)oip{x) 



(11) 



where Sf is the action for free fermions with momentum dependent mass, 



V 



dt dx dyip(x, t) [iS(x — y) r y fi d^ — m[x — y)] ip(y, t) . 



(12) 



In the subsequent calculation we will treat e<p as small and expand to second order 
the logarithm in the exponential of the norm produced by integrating over ip and ip. This 
results in a series of diagrams. The leading order is independent of and contributes a 
constant factor to the norm which will cancel. Represented by Fig. 1, the term quadratic 
in <p contains a fermion loop and is written, 

o2 



e 
2~ 



dx dy <p{x)<i){y)TY [y°D + (x - y)-y°D_(y - x)] , 
where the equal time propagators are, 



(13) 



D±(k) 



o i_ V + m(k) 



7 U ± 



E(k) 



(14) 



and E 2 (k) = k 2 + m 2 (k). We will see later that keeping only second order terms in e<p is 
sufficient to locate the best variational state. The contribution to the norm from fermion 
integration is, 
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(m| exp \^—ie J dx(f)(x)ip^(x)il ! (x)J \m) = 

const, x exp J dx dy d(p(x) X(x — y) d(f>(y) + O ((e0) 3 )^ , (15) 

where, 

x(fc) - * y * i 1 ^(p-fc) — J • (l6) 

The norm then can be written as Z = ZaZ^Z^, where Za and Z^ are constants which cancel 
in all calculations and 

Z^ = J D<P exp j dxdyd^x) S{x - y) d<f>{y)j . (17) 

The "effective action" for <fi is given by 

e 2 

S(k) = G- l {k) + -X{k) . (18) 

71 

We now compute the expectation value of the Hamiltonian in the variational state. Min- 
imization of this quantity with respect to the variational parameters will yield the ground 
state of the theory. In doing so one should be careful with properly regularizing the ultra- 
violet singular operators that enter the Hamiltonian, 

H{x) = ^E 2 (x) + #(x)7i (d + ieA(x)) i/;(x) . (19) 

The term to keep an eye on is the interaction term. In particular the current j(x) = 

^(x) r yijj(x) has to be regularized using the point-splitting technique [^]. This point-splitting 

generates an additional term in the Hamiltonian which is quadratic in A |7]|1 . 

A 2 - e 2 ie 

ii>Drfi> = -— + % : V^yV : -e A : V>7 V : -^-A 2 + —OA , (20) 

Z7T Z7T ATT 

where A is the UV cutoff. The first term is a constant and it does not affect our calculation. 
The last term is a total derivative and for this reason will also be dropped. We are therefore 
concerned with the calculation of the expectation of the regularized Hamiltonian, 

H = j dx f^E 2 (x) + : i$(x)'yidip(x) : +eA(x) : i>(x)^(x) : +^-A 2 (x)\ ■ ( 21 ) 



We proceed to calculate the expectation of ([H]) in the manner outlined above for the calcu- 
lation of the norm. We find, 



(i dx : ip(x) n fidip(x) :) = lim 



dk 



k 7 



A 2 



A ^°° \J A 27r ^k 2 + m 2 (k) 2tt^ 
2 ' dk 



(e I dz$(z)<yAil>(x)) = —J ^S~\k) , 

(^ldxA\x)) = ^j^[G(k) + S-\k)]. 



(22) 

(23) 

(24) 
(25) 



Functional variation of the expectation of the Hamiltonian with respect to the parameters 
m(p) and G~ 1 (p) yields two minimization equations, 



5m(p) 
5(H) 



d ±. (c-Hk) + -) S- 2 (k) 5X{ - k) - p2m ^ P) 
8tt V vr; >5m{p) (p 2 + m?{p)f /2 



, 



s-\p) 



1 - -G 2 (p)) X 2 {p) - — X(p)G(p) 
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7T 



. 



(26) 
(27) 



5G-\p) 4 

Since 5X/5m oc m, the first equation is solved for m(p) = 0. The second minimization 
equation has the solution 



-G(p) + X-\p)) = e l + X- 2 {p) , 

7T / 7T 



(28) 



l/|p|, and we can write the 



where X(p) is defined in (|i~6|) . With m(p) =0 we have X(p) 
solution as, 

2 \ ^ 2 

-G{p) + \p\) =-+ P 2 . 

71 J 71 

Summarizing, the result of our variational calculation is the vacuum wavefunctional 
the parameters 
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m(p) = 0, G(p) = -\ J-+ p 2-\p\ . 
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(29) 



with 



(30) 



The calculations in this section rely upon the perturbative expansion of the "effective action" 
for in (|T5D. To improve on this we would have to resum this expansion. However, for 



S 



m(q) = such a re-summation is trivial. It can be shown explicitly that all diagrams higher 
than second order vanish and therefore our result is exact at least in the sense that the 
solution we found is at least a local minimum. We will see in the next section that the result 
is indeed exact. 



III. BOSONIZED FORMULATION OF THE SCHWINGER MODEL. 

As is well known the Schwinger model is exactly solvable by the bosonisation technique 
||. In this section we will discuss the variational calculation in this bosonized form. We 
shall show that the solution obtained by the variational method in the bosonized form is 
exact. Further, we shall show that this solution is equivalent to the wavefunctional obtained 
in the previous section constructed in the fermionic Hilbert space. 

The bosonized form of the Hamiltonian density in the temporal gauge (A = 0) reads, 



H(x)= 1 -E\x) + Up(x)-j=A(x: 



1 2 

2 l d x(x)} 



which is supplemented by a Gauss' law constraint on physical states, 



(31) 



Gb{x) = E(x) = x(x 







7T 



(32) 



where E(x) and A(x) = A\(x) are the electric and photon fields as defined in the previous 
section, xi x ) is the boson field, and p(x) its canonical momentum. Since the Hamiltonian 
is quadratic in terms of bosonic fields and their conjugate momenta, the exact ground state 
wave functional can be readily found, 



*o[Ax] =exp 



J dxx{x)A(x) - - I dxdyx{x)^ l {x-y)x{y) 



(33) 



with Ti~ 1 (k) = \/k 2 + e 2 /ir and the suffix indicates the exact ground state. 

We now want to show that this result is recovered within the variational calculation. We 
base our variational Ansatz upon a Gaussian, 



cxp 



dxdyA(x)G 1 (x - y)A(y) - ^ / dxdyx{x)T, x {x - u)\ f //) 



(34) 
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with G{x — y) and E as variational parameters. Projecting this onto the gauge invariant 
subspace with the help of the Gauss' law eq.(^) we have the gauge invariant trial state, 

^b[A, X } = J D4> exp -~ j dxdyA^(x)G- 1 (x-y)A'f'(y) 

~ j dx dy x{x)t' 1 (x - y)x(y) + ^= f dx d(f>(x)x(x) ■ (35) 



Proceeding with the variational calculation in the usual way, we find the expectation 
value of the Hamiltonian to be, 

-G-\k)t(k) + ( tifyGik))" 1 + — + k 2 t(k)G~\k) + ^-G(k)t(k) 

T V / 7T IT 2 



lJM(y' {k) + e + *l)Y-i {k) , (36) 



where 
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Y(k) = E-\k) + —G(k) . (37) 

7T 

Variational minimization of (|36| ) with respect to Y(k) yields the required solution, 

Y 2 (k)= (£-\k) + ^-G(k^j =(k 2 + ^j. (38) 

After substituting into (|35|) and integrating over <fi, we find that the best variational state 
is indeed (|33|), so the variational calculation gives the exact result. 

Note that our initial parameterization of the trial state turned out to be somewhat 
redundant. The energy minimization determines only Y but not G and S separately. It is 
in fact true that integrating over <ft in equation (|3^) we obtain a Gaussian in x with width Y. 
This redundancy however can be used to establish the equivalence of this calculation with 
the calculation in terms of the fermionic representation of the previous section. The point 
is that we can choose £ to be the width of the ground state in the free massless bosonic 
theory, 

= \fk? . (39) 
10 



The exact ground state is then represented as a gauge projected product of the vacuum wave 
functional of a theory of a massless boson and a Gaussian functional of the gauge potential 
A. In fact for this choice of S, eq.(|38|) gives the same result for the width of the gauge field 
Gaussian G as the fermionic solution eq. fl30|) . Also the free massless fermion by bosonisation 
is equivalent to the free massless boson. Therefore the fermionic part of the variational state 
of the previous section for m = is equivalent to the free massless state for the bosonic field 
X- We therefore see that our variational solution of the previous section is equivalent to the 
solution obtained here, and is therefore an exact ground state. 

IV. CONCLUSIONS 

In this paper we have presented a variational analysis of the massless Schwinger model. 
We have shown that there exists a simple and natural way to extend the gauge invariant 
variational Ansatz of |I| to systems with fermions. We have also found that in this simple 
and well understood system our variational Ansatz in fact includes the exact ground state 
and therefore reproduces known exact results. 

A noteworthy feature of the best variational state is that the dynamical mass in the 
fermionic sector is not generated, so that in terms of fermionic part the state seems to be 
trivial. In fact in 2 dimensional theory this is not entirely surprising. If we were to consider 
a theory with more than one flavour, this would have been an immediate corollary of the 
Coleman theorem. Since a continuous (chiral) symmetry can not be broken in 2 D, the state 
with nonvanishing dynamical mass should be strongly energetically disfavoured. In the one 
flavour theory which we are considering here, such an a priori argument does not apply since 
the axial U(l) symmetry is anomalous. Nevertheless it is not unnatural that in terms of the 
dynamical mass the ground states in one flavour and multiflavour theories are similar. 

On the other hand in the one flavour case the fermionic bilinear condensate ipip should 
be nonvanishing precisely due to the same anomaly. It is therefore interesting to see whether 
our best trial state leads to such a nonvanishing condensate. 
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It is simplest to calculate the condensate first in the bosonized version of the theory. The 
bosonization rules prescribe the identification M, 



ip{x)^{x) 



ce 



7T 



(40) 



where M e i^ means normal ordering with respect to the free boson field of mass ej '-y/7T, and 
the prefactor is c = exp (j)/2 , k, with 7 being Euler's constant We need the average of 
(fiO|) over the wavefunctional fl35|) . After integrating over A and x we obtain 

1 



ce 



= N e /^ / exp 



( dydzdcf>(y) \G- l {y-z) + -H{y-z)) d<f>(z) 



+e / dyt{x-y)d<l>{y)-TrE{0) 



(41) 



Using the definition a further integration over yields 



(ip(x)ip(x)) 



ce 



Ne/y/Z ex P 1 ( x - &)] 



(42) 



The exponential is UV divergent, but it cancels exactly with the normal ordering factor, so 
we finally obtain 

e exp (7) 



(ip(x)ip(x)) 



2tt 3 / 2 



(43) 



which is the exact result 11 



The same calculation can be performed in the fermionic formalism. We need to calculate 

1 



Dtjj Di/jijj(x)ijj(x) exp 



dy dz d(j){y)G (y — z)d<f)(z 



+i / dt dx ip(x, t)7 M (idfj, — ea^) tp(y, t) 



(44) 



with a (x, t) = 4>(x)5(t) and a\(x, t) = 0. Here we have used the path integral representation 
of the average over the massless fermion vacuum state eq . (|i~l~l) . To perform the integration 



over the fermions we can use the results of |12| . The fermionic integral is non- vanishing 
only for background fields with instanton number ±1. Decomposing a M = 8^6 — e^d^ip, we 
obtain 
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(ip(x)ip(x)} = const, x I Dcj) exp 
1 



dy dz d(f>(y) S(y — z) d<f>(z) 



2ty 



dy \og(x - y) d<p(y) 



(45) 



Here we have used 



0<p(x, t) = -d(j)(x) S(t) <f(x, t = 0) = — [ dy log{x - y)d<j>(y) . 

2ir J 



(46) 



With S given by eq.(J5S|) this coincides with eq.(^TJ) and so the subsequent integration 
over (ft again reproduces eq. (|4"3"D . 

We have thus seen that in this simple model the gauge invariant variational Gaussian 
approximation works well. Extension of this approach to four dimensional QCD with chiral 
symmetry breaking is the next challenging task. 
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